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A crossover between logarithmic and exponential temper- 
ature dependence of the conductance (weak and strong lo- 
calization) has been observed in ultrathin films of metals de- 
posited onto substrates held at liquid helium temperatures. 
The resistance at the crossover is well denned by the onset 
of a nearly linear dependence of conductance on thickness at 
fixed temperature in a sequence of in situ evaporated films. 
The results of a finite size scaling analysis treating thickness 
as a control parameter suggest the existence of a T = quan- 
tum critical point which we suggest is a charge, or electron 
glass melting transition. 

PACS numbers: 71.30.+h, 05.70. Jk, 72.15.Rn, 73.50.-h 



Recently there has been considerable interest in electri- 
cal transport in two dimensional (2D) conductors. Since 
the work of the so-called "gang of four" it has been the 
conventional wisdom that in two dimensions weak disor- 
der will localize electrons making it impossible for true 
metallic behavior to occur M. Although the classical 
work focused on the problem of noninteracting electrons, 
it has been widely accepted that its conclusions were ap- 
plicable to interacting systems as well. This view has re- 
cently been strongly challenged by experimental studies 
of ultrahigh mobility low carrier density Si MOSFETs 
[^]-||] which have revealed in an elegant way a metal- 
insulator transition (MIT) in a two dimensional system. 
Dobrosavljevic et al. |6| in light of the work of Refs. (^-^] 
recently reconsidered the scaling theory of the MIT in two 
dimensional interacting systems and asserted that the ex- 
istence of a 2D metal-insulator transition does not contra- 
dict any general principles. Another recent development 
in the physics of two dimensional conductors was the ob- 
servation of Hsu et al. of a change in the sign of the 
magnetoconductance of ultrathin films with strong spin- 
orbit interactions. This was correlated with a change 
in the temperature dependence of the conductance from 
logarithmic to exponential with decreasing temperature, 
i. e. the crossover between weak and strong localization, 
an effect which had been observed previously MM . 

Here we present the result of a detailed study of the 
dependence of the conductance on thickness (d) and tem- 
perature in ultrathin Pd j^] and Bi films, metals with 
strong spin-orbit interactions. The success of a finite-size 
scaling analysis of these data leads us to suggest that 
there is a quantum critical point (QCP) in these films 
with the Drudc conductance of the film as the control 



parameter. A QCP is found at zero temperature when a 
ground state is controlled by an external parameter [fl0[ . 
It is our suggestion that the QCP we have found may 
be a charge or electron glass melting transition. This is 
motivated by the possibility that features of the descrip- 
tion of vortex hopping in a two dimensional disordered 
superconductor may be similar to electron hopping in a 
two dimensional disordered conductor . Although the 
vortex-glass phase transition in two dimensions is sup- 
pressed to T = 0, a 2D vortex glass correlation length 
which diverges in the limit of zero temperature has been 
reported pq |. This implies that there may be a QCP in 
this system as well. Finally there is ample evidence that 
high resistance films of metals whose conductance is gov- 
erned by hopping exhibit glass-like properties, whereas 
low resistance films whose conductance is logarithmic in 
temperature, and which are not superconducting are not 
glass-like fjj. 

Our application of finite-size scaling to the weak-to- 
strong localization crossover was motivated by the obser- 
vation that although the crossover appears gradual as a 
function of temperature as in Ref. 0, there is actually 
a very well-defined crossover conductance, G cr , which is 
particularly sharp especially when the thickness, d, is the 
independent variable. We find G cr ~ (26 ± 3fcf2) _1 and 
G cr ~ (56 ± 6fc£l)~ 1 for Bi and Pd films respectively. 
Although G cr is different in these materials, the criti- 
cal exponents of the scaling analysis are nearly identical 
within experimental uncertainty. 

Film thicknesses ranged from 8 to 20A. Bismuth films 
were grown on a 9A pre-deposited layer of amorphous Ge 
(a-Ge), whereas Pd films were grown directly on glazed 
ceramic substrates ||. Depositions were carried out in 
situ under UHV conditions (< 10~ 9 Torr) with substrate 
temperatures held below 15K. The ultra-high vacuum 
environment of the samples was sustained over an ex- 
tended period so that film thickness could be built up 
by a scries of sequential depositions without contamina- 
tion. For a given set of films the deposition conditions 
(evaporation rate, source and substrate temperature, and 
background pressure) were fixed to avoid changes in mor- 
phology from film-to-film in a sequence. Resistance was 
measured between deposition increments using standard 
four-probe techniques with films current-biased at values 
less than InA. Current- voltage characteristics were linear 
up to currents at least ten times this value. 

It is generally accepted that films of metals deposited 
on a-Ge substrates held at helium temperatures are con- 
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tinuous and homogeneous for two reasons: first, measur- 
able conductance is often found at the order of mono- 
layer coverages. Second, at very early stages of growth, 
the conductance becomes a linear function of thickness. 
Palladium films are known to wet glass substrates, also 
yielding films with measurable conductances at mono- 
layer coverages. Close examination of our data indicates 
that a linear dependence of the conductance on thickness 
is never really found for films at the very earliest stages 
of their growth, even somewhat beyond monolayer cov- 
erage. An analysis of the onset of nearly linear behavior 
with increasing thickness is the subject of this work. 

Figure |l| shows the conductance, G, of a Bi film as a 
function of d at several representative temperatures, and 
the derivative of G(d) with respect to d at 3.3K. These 
data are obtained by sorting data on G(T) obtained from 
films of different thicknesses, and considering them as a 
function of d at fixed temperature. The derivative is then 
computed numerically. Graphs of the derivative, of which 
only one example is shown, exhibit a pronounced change 
at a temperature dependent thickness d c which always 
occurs at a constant conductance G cr . Films with thick- 
nesses greater than d c exhibit a nearly linear dependence 
on thickness. The values of the critical conductance for 
Bi and Pd films, G cr ~ (26±3fcfi)" 1 , and (56±6fcO)-\ 
respectively, are very close to e 2 /h and e 2 /2h. 

Explanations of G cr being a consequence of a struc- 
tural change occurring at a particular thickness can be 
ruled out. If such were the case, the boundary between 
the linear and non-linear dependences would occur at a 
fixed thickness independent of temperature rather than 
at a fixed conductance G cr . It should also be mentioned 
that the conductivity near treshold is not described by a 
percolation model with the conductance given by a power 
law in {d — dp), where d p is the thickness at the threshold 
for conductance. In this approach the average thickness 
is assumed to be proportional to the surface coverage 
JTij . It is actually given by an exponential in l/d a , al- 
though the power a is somewhat uncertain because of 
the limited range of the data. A value of a = 1/2 would 
be consistent with variable range hopping |l5| . 

Further support of the view that the crossover is an in- 
trinsic effect comes from the temperature dependence of 
the conductance at fixed thickness, which is not shown. 
At low temperatures, and in the thinnest films at all 
temperatures, the conductance is activated. For resis- 
tances below 50 — lOOfcri it is given by the well-known 
functional form of the soft Coulomb gap model G oc 
exp[-(T /T) 1 / 2 ] Here T is related to the local- 

ization length, Ll oc , through fcsTo = e 2 /kLl OC) where 
k is the dielectric constant. For resistances above 50- 
100kf2 the exponent is closer to 4/5 rather than 1 /2 which 
may be evidence of collective variable range hopping of 
charges, or may be the result of a crossover from variable 
range to fixed range hopping with increasing resistance. 
At high temperatures for all but the most resistive films 
the conductance varies linearly with In T consistent with 
the expectations of the electron interaction picture p6[. 



Because the crossover in temperature between logarith- 
mic and activated conductance is not as sharp as the 
thickness-dependent crossover, it has not been used to 
determine G cr . The fact that the crossover is seen at the 
same G cr in both the temperature and thickness depen- 
dence of the conductance suggests that the effect is not 
some manifestation of percolation on an atomic scale fHj ] 
. We have also observed a change in the sign of the mag- 
netoresistance in Bi films from negative to positive with 
increasing thickness which is qualitatively consistent with 
Ref. [Q. However, it does not occur at the same value 
of the resistance as the crossovers described above. This 
delicate issue is currently under study. 

The above observations led us to investigate the possi- 
bility of a quantum critical point (QCP). In this scenario 
for T > there would be a line of crossovers at d = d c (T), 
terminating in a critical point at T — 0. It is sufficient 
for a QCP that the finite temperature boundary be a 
crossover. It is not necessary that it be a line of first or- 
der transitions ]l8| ] . Because we can determine G cr from 
G(d) with precision we take thickness d as the control pa- 
rameter, and define d c as a temperature dependent crit- 
ical thickness. The use of a temperature-dependent pa- 
rameter will be seen to be a convenience, with the Drude 
conductance at high temperatures, which is temperature 
independent, being the actual control parameter. 

The distance from the critical thickness, 6 = d — d c 
is then the relevant field. Standard arguments ||io| , [l8| , [l9f 
suggest that the scaling form for a two dimensional sys- 
tem is 

G{6,T) = G cr F{ST- 1 !™) (1) 

Here F(x) is a universal scaling function such that 
F(0) = 1, v is the coherence length exponent, and z 
is the dynamical critical exponent. To proceed with the 
analysis we rewrite Eq. |l| 

G(S,t) = G cr F(St) (2) 

and define t = T~ x l vz . The parameter t(T) is treated as 
an unknown variable to be determined at each tempera- 
ture to obtain the best collapse of all the data p9| . The 
exponent vz is then found from the temperature depen- 
dence of t, which must be a power law in temperature 
for the procedure to make physical sense. We should 
note that the present scaling procedure does not require 
detailed knowledge of the functional form of the temper- 
ature dependence of the conductivity, or prior knowledge 
of the critical exponents. It is simply based on the data 
which includes an independent determination of d cr at 
each temperature. 

The results of the analysis for both Bi and Pd are 
shown together in Fig. ||. The crossover thickness d c (T), 
determined from the data like that of Fig. [j], increases 
logarithmically with decreasing temperature, as shown 
on Fig. |^. The temperature dependence of the parame- 
ter t, shown on Fig. ^is consistent with a power law, with 
vz = 6.9 ± 0.7 and 7.2 ± 1.0 for Bi and Pd, respectively. 
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Although our analysis is carried out for convenience 
with an apparently temperature-dependent critical pa- 
rameter, it can be cast in terms of a temperature in- 
dependent one. Adding metal sequentially to a quench- 
condensed film controls disorder through its smoothing of 
the random potential. What probably also occurs is an 
increase in the carrier concentration, and thus an increase 
in the screening. These effects are measured directly by 
the conductance, and indirectly by the thickness. Al- 
though G is always a temperature dependent quantity 
its value at high temperatures is close to the Drude con- 
ductance Gb i which is proportional to the thickness and 
constant for a given film. A natural control parameter 
for the QCP would be Go, with the relevant field being 
5d = \Gd — G cr | . Because all the films we used in the scal- 
ing analysis exhibit InT behavior in their conductances, 
with G oc d at high temperatures, the control parameter 
used in our analysis, 5, is effectively the same as 5d pro- 
vided that d c (T) is a logarithmic function of temperature 
as is demonstrated in Fig||. In other words, temperature 
dependence of d c is not a consequence of an underlying 
finite temperature phase transition, it only reflects the 
fact that d c is related to G cr through a temperature de- 
pendent factor. 

The crossover will always occur in films which arc effec- 
tively two dimensional at nonzero temperature because 
the inelastic scattering length, li n , will always be greater 
than the film thickness. One can crudely estimate £;„ by 
equating it to the localization length at the crossover con- 
ductance |^(| ■ This is tantamount to matching the length 
scales of weak and strong localization at the crossover. 
Taking the dielectric constant to be that of Ge, k = 16, 
one finds U n ~ 100A in the thinnest films at the crossover. 
The inelastic scattering length, ij n , is a power law in 
temperature proportional to T~ x / 2 which would increase 
from ~ f 00A with decreasing temperature. Since d c in- 
creases as oc InT from ~ 12 A over the measured temper- 
atures, one would expect d c to remain smaller than the 
localization length at all nonzero temperatures. Exami- 
nation of the data of Fig. 2 suggests that the temperature 
dependence of <i c may be saturating at the low end of the 
temperature range, which may mean that <i c does not in- 
crease indefinitely and is bounded at temperatures below 
the measured range. 

In previous work done by our group Q and others , 
the data for the thickness and temperature dependences 
of sequences of films have been scaled by adjusting only 
the temperature axis. This procedure did not include a 
specific recognition of the crossover phenomena shown in 
Fig. [IJ and involved forcing a functional form on high 
temperature data which was slightly different from that 
of localization theory 0] . 

The QCP separates films whose conductances exhibit 
logarithmic from those with exponential dependence on 
temperature at nonzero temperatures. Although both 
types of films are nominally insulating, the analysis would 
suggest that they have either fundamentally different 
ground states or perhaps different dynamical behavior. 



The existence of a serious distinction between insulators 
exhibiting logarithmic and exponential dependence on 
temperature was first noted in the work of Burns and 
Chaikin |H| on thermopower in Pd and Pd-Au films. 
They found a sharp crossover in the thermopower from 
metallic-like to insulating-like behavior at a resistance 
near h/e 2 . There is another important difference be- 
tween films whose conductances exhibit a logarithmic de- 
pendence on temperature and those with an exponential 
dependence. The former appear never to be glass-like 
in their behavior, whereas the latter are unequivocally 
glass-like [^3| . The evidence for this is the observation of 
long relaxation times and memory effects in a capacitive 
charging experiment. Glass- like behavior disappeared 
above a characteristic temperature roughly correspond- 
ing to the crossover from hopping to a logarithmic de- 
pendence of the conductance on temperature. 

The magnitude of the critical exponent products we 
have found (vz ~ 7) seem large. In the absence of a 
detailed theory, it is hard to comment on this. However 
the value is not very different from the exponent prod- 
ucts vz ranging between 4 and 10 found by Dekker et 
al. jj3 in their study of the vortex-glass phase transi- 
tion suppressed to zero temperature in a 2D supercon- 
ductor. Thus, the large values of critical exponents in 
the present work would be plausible if charges in disor- 
dered ultrathin films and vortices in two dimensional su- 
perconductors behaved similarly. This line of argument 
leads us to suggest that the QCP inferred from the anal- 
ysis we have presented may be the melting of a charge or 
an electron glass, with the higher-resistance phase being 
glass- like with activated transport, and the lower resis- 
tance phase being a liquid-like with logarithmic correc- 
tions to the classical conductivity from electron-electron 
interaction effects. 

We gratefully acknowledge numerous useful discussions 
with Myron Strongin, Daniel Grupp, Boris Shklovskii, 
Anatoly Larkin and Subir Sachdev. This work was sup- 
ported in part by the National Science Foundation under 
Grant No. NSF/DMR-9623477. 



[1] E. Abrahams, P. W. Anderson, D. C. Licciardello, and 

T. V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 
[2] S. V. Kravchenko, W. E. Mason, G. E. Bowker, J. E. 

Furneaux, V. M. Pudalov, and M. D'lorio, Phys. Rev. B 

51, 7038 (1995). 
[3] S. V. Kravchenko, G. V. Kravchenko, J. E. Furneaux, 

V. M. Pudalov, and M. D'lorio, Phys. Rev. B 50, 8039 

(1994). 

[4] S. V. Kravchenko, D. Simonian, M. P. Sarachik, Whitney 
Mason, and J. E. Furneaux, Phys. Rev. Lett. 77, 4938 
(1996). 

[5] Dragana Popovic, A. B. Fowler, and S. Washburn, Phys. 



3 



Rev. Lett. 79, 1543 (1997). 
[6] V. Dobrosavljevic, Elihu Abrahams, E. Miranda, and 

Sudip Chakravarty, Phys. Rev. Lett. 79, 455 (1997). 
[7] S.-Y. Hsu and J. M. Valles Jr., Phys. Rev. Lett. 74, 2331 

(1995) . 

[8] Y. Liu, B. Nease, K. A. McGreer and A. M. Goldman, 

Europhys. Lett. 19, 409 (1992). 
[9] S. J. Lee, J. B. Ketterson, and Nandini Trivedi, Phys. 

Rev. B 46, 12695 (1992). 
[10] S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Shahar, 

Rev. Mod. Phys. 69, 315 (1997). 
[11] M. P. A. Fisher, T. A. Tokuyasu, and A. P. Young, Phys. 

Rev. Lett. 66, 2931 (1991). 
[12] C. Dekker, P. J. M. Woltgens, R. H. Koch, B. W. Hussey, 

and A. Gupta, Phys. Rev. Lett. 69, 2717 (1992). 
[13] G. Martinez-Arizala, C. Christiansen, D. E. Grupp, N. 

Markovic, A. M. Mack, and A. M. Goldman, Phys. Rev. 

B57, R670 (1998). 
[14] D. Stauffer, Phys. Rep. 54, 1 (1979). 
[15] A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1975). 
[16] B. L. Altschuler, A. G. Aronov and P. A. Lee, Phys. Rev. 

Lett. 44, 1288 (1980). 
[17] Zhoa-Quing Zhang and Ping Sheng, Phys. Rev. B 44, 

3304 (1991); J. C. Phillips, Phys. Rev. B 42 2825 (1990); 

Y. Mcir, A. Aharony, and A. B. Harris, Europhys. Lett. 

10, 275 (1989). 

[18] Elihu Abrahams and Gabriel Kotliar, Science 274, 1853 

(1996) . 

[19] Daniel E. Grupp and Allen M. Goldman, Phys. Rev. Lett. 

78, 3511 (1997). 
[20] Patrick A. Lee and T. V. Ramakrishnan, Rev. Mod. 

Phys. 57, 287 (1985). 
[21] M. J. Burns and P. M. Chaikin, Phys. Rev. B.10, 5924 

(1983). 

FIG. 1. G vs. d for Bi films at 0.7K (o), 1.5K (o), 3.3K (•), 
7.0K (A) and 15. OK (v). A fit to the linear part of the 3.3K 
curve is shown. Also plotted is the derivative of the 3.3K 
curve (x). The horizontal dashed line represents the slope 
of the linear fit. A sharp bend in the derivative is observed 
very close to the thickness at which G(d) = e 2 /h. Efforts to 
fit these curves to a percolation model, G cx (d — dp) 1 were 
unsuccessful. 

FIG. 2. Bi (o) and Pd (•) data collapsed together. The 
poor collapse of the Pd data in the lower right is from the 
highest temperature curves (T > 1QK), indicating that the 
scaling breaks down far from the T — transition. 

FIG. 3. d c vs. T for Bi (o) and Pd (•). d c increases slowly 
(logarithmically) with decreasing temperature. 

FIG. 4. T vs. t(T) for Bi (o) and Pd (•). The slopes of the 
power law fits yield vz = 7.2 ± 0.5 and vz = 6.9 ± 0.5 for Bi 
and Pd, respectively. 
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